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•  Based  on  the  field  measurements  and  the  modeling  effort  performed  ,  it  is  clear  that  the  tidal  cycle,  the 
wind  setup,  and  the  backwater  effects  dominate  the  hydrologic  response  of  coastal  low-gradient 
watersheds  to  rain  storms.  Therefore,  in  order  for  any  predictive  tool  to  produce  accurate  and  reliable 
predictions  it  has  to  incorporate  these  factors. 

•  A  well-trained  Artificial  Neural  Network  and  Non-Parametric  Regression  models  have  been  developed 
and  validated.  The  models  were  able  to  reproduce  the  complex  and  non-unique  stage-discharge 
relationship  observed  at  the  outlet  of  tidal  low-gradient  watershed.  It  should  be  emphasized  that 
standard  statistical  tools  completely  failed  to  capture  this  complex  stage-discharge  pattern. 

•  The  comparison  between  physically-based  and  lumped  hydrologic  models  applied  to  the  Goodwin  Creek 
watershed  showed: 

o  The  predictive  ability  of  both  modeling  approaches  (when  well  calibrated)  for  mid-size 
catchment  was  surprisingly  similar. 

o  Their  response  to  variations  in  the  temporal  sampling  of  rainfall  was  also  similar. 

o  Their  response  to  variations  in  the  spatial  density  of  rainfall  information  is  currently  being 
investigated. 

•  An  extensive  modeling  effort  was  prefonned  to  establish  guidelines  on  the  temporal  and  spatial  rainfall 
sampling  requirements  for  hydrologic  prediction.  A  resolution  frequency  of  up  to  one  hour  seems 
adequate  for  practical  purposes  of  predicting  peak  runoff  discharges. 

•  An  equipment  grant  from  the  State  of  Louisiana  was  received  in  support  of  the  data  collection  and 
modeling  efforts  of  this  project.  The  equipment  purchased  included  telemetric  units  of  precipitation, 
stage,  runoff  discharge,  total  weather  stations,  and  rainfall  distrometer. 

•  A  supplement  funding  through  the  Louisiana  Transportation  Research  Center  was  also  received  to 
develop  a  comprehensive  drainage  plan  for  the  low-gradient  Vermilion  watershed. 

•  A  proposal  to  the  State  of  Louisiana  Board  of  Regents  was  recently  funded  to  set  up  an  experimental 
study  to  investigate  the  uncertainty  of  radar  rainfall  estimates  and  effects  on  hydrologic  flood  prediction. 

(5)  Technology  transfer. 

•  Work  closely  with  the  Coastal  Hydraulics  Laboratory  (CHL)  of  the  Engineering  Research  and 
Development  Center  (ERDC)  hydrologic  research  group.  The  research  activities  of  this  project  include 
continuous  improvements  and  development  to  the  GSSHA  model. 

•  The  development  work  for  GSSHA  is  coordinated  with  the  University  of  Connecticut  and  BYU. 
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c.2)  List  of  Illustrations 

Figure  1:  Layout  of  Isaac  Verot  catchment  along  with  the  equipment  installed.  The 
rectangular  grid  represents  the  4x4  km2  Stage  III  radar  pixels. 

Figure  2:  An  example  of  GIS  analysis:  Soil  classifications  maps  within  the  Vermillion 
Watershed. 

Figure  3:  Relationship  between  stage  and  discharge  in  flat  watersheds  subjected  to  tidal 
influence  is  complex  and  quiet  challenging.  Data-driven  models,  such  neural  networks, 
can  reproduce  these  relationships  with  relatively  high  accuracy 

Figure  4:  Observed  and  simulated  hydrograph  for  the  three  models  during  calibration 
period 

Figure  5:  Observed  and  simulated  hydrographs  for  the  three  models  during  a  validation 
period 

Figure  6:  Error  in  rainfall  and  runoff  due  to  rainfall  temporal  resolution. 

Figure  7:  Error  in  rainfall  and  runoff  due  to  changes  in  spatial  resolution  of  rainfall 

Figure  8:  The  2 1-km  Goodwin  Creek  Experimental  Watershed  in  northern  Mississippi 
is  selected  for  radar-rainfall  runoff  analysis  because  of  its  dense  coverage  with 
precipitation  and  surface  runoff  measuring  stations. 

Figure  9:  MPE  radar  estimates  were  found  to  have  significant  bias  during  most  of  the  1 1 
storms  selected  for  this  study.  Storm  bias  factors  for  each  pixel  covering  the  watershed 
were  computed  as  the  ratio  of  total  radar  to  total  gauge  accumulations.  Bias  factors 
showed  significant  variability  from  one  pixel  to  another. 
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c.3)  Statement  of  the  Problem  Studied 

The  goal  of  the  project  is  to  assess  the  ability  to  make  hydrologic  predictions  for  low- 
gradient  watersheds,  and  to  determine  the  challenges  that  such  watersheds  pose  to 
hydrologic  models.  The  Vermilion  Watershed  and  its  sub-catchments  in  southwest 
Louisiana  were  selected  to  test  widely  used  numerical  hydrologic  models.  The  Vermilion 
watershed  has  a  very  low-gradient  topography  and  very  complex  fluvial  systems  and  tidal 
channels.  Making  accurate  hydrologic  predictions  for  these  systems  is  a  challenging  task 
because  of  the  dominant  role  of  hydraulics. 

Research  activities  of  this  project  focused  on  the  following  research  problems: 

i.  Collection  and  analysis  of  hydrological  and  meteorological  measurements  in  a 
low-gradient  watershed  in  southwestern  Louisiana. 

ii.  Development  of  a  comprehensive  GIS  database  for  the  watershed,  including  the 
topographical  and  the  hydrological  characteristics;  e.g.  land  use,  soil  types, 
storage,  infiltration  parameters  and  roughness. 

iii.  Development  of  data-driven  models  (Artificial  Neural  Networks  and  non- 
parametric  local  regression)  for  predictions  of  complex  discharge-stage  relations 
in  low-gradient  streams. 

iv.  Performance  evaluation  of  GSSHA  modeling  system  compared  to  other  widely 
used  hydrologic  models  such  as  MIKE-SHE  and  HEC-HMS. 

v.  Assessment  of  sensitivity  of  runoff  predictions  to  spatial  and  temporal  rainfall 
sampling. 

vi.  Using  radar  rainfall  information  for  hydrologic  modeling  with  GSSHA. 

c.4)  Summary  of  the  Most  Important  Results 

i.  Collection  and  analysis  of  hydrological  and  meteorological  measurements: 
Through  preliminary  applications  of  CASC2D,  it  became  apparent  that  the  existing 
network  of  monitoring  stations  for  rainfall,  stage,  and  discharge  was  not  adequate  to 
setup  and  calibrate  a  numerical  hydrologic  model.  The  investigators  were  able  to 
purchase  equipment  to  collect  such  data  though  state  and  local  funding.  The 
investigators  determined  that  it  would  be  more  feasible  to  focus  on  a  mid-size  sub¬ 
catchment  of  the  Vermilion  Watershed,  namely  Isaac  Verot.  The  instrumentations 
were  placed  in  this  sub-catchment  with  an  adequate  density  such  that  thorough 
hydrologic  analysis  can  be  conducted.  A  brief  summary  of  the  acquired  equipment 
is  listed  below: 

•  13  dual-rain  gauge  stations  with  automatic  data  loggers. 

•  A  total  of  five  discharge  gauges  installed  at  the  watershed  outlet  and  at  interior 
sites.  All  gauges  are  YSI  Sontec  Side  and  bottom-looking  acoustic  devices. 

•  One  weather  station  that  measures  solar  radiation,  wind  speed  and  direction,  soil 
moisture,  relative  humidity,  air  temperature,  and  barometric  pressure. 

•  One  Joss-Wldvogel  acoustic  rainfall  distrometer. 

•  All  monitoring  stations  are  equipped  with  telemetric  units. 
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4  Km 


O  Dual-rain  gauge  site 
U  J-W  Disdrometer  and  dual-rain  gauge  site 
A  Weather  station  and  dual-rain  gauge  site 
^  Discharge  gauge  site 
O  USGS  rain  gauge 

Figure  1:  Layout  of  Isaac  Verot  catchment  (southwestern  Louisiana)  along  with  the 
equipment  installed.  The  rectangular  grid  represents  the  4x4  km2  Stage  III  radar  pixels. 
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ii.  Development  of  a  comprehensive  GIS  database: 

A  comprehensive  GIS  database  for  the  watershed  was  developed.  For  accurate  and 
consistent  representation  of  the  watershed  drainage  area,  the  US  Geological  Survey 
(USGS)  hydrologic  unit  boundaries  were  adopted.  The  environmental  analysis 
system  BASINS  developed  by  the  Environmental  Protection  Agency  (EPA)  was 
used  to  provide  these  boundaries  in  a  GIS  geo-referenced  format.  BASINS 
database  was  also  integrated  in  the  GIS  database  for  land  use  information,  locations 
of  existing  meteorological  and  flow  measurements  stations,  and  detailed  stream 
network  within  the  watershed.  The  topography  of  the  watershed  was  determined 
from  the  30-meter  USGS  Digital  Elevation  Models  (DEM).  Soil  types  and 
characteristics  were  compiled  and  incorporated  into  the  database.  In  addition, 
recent  high-resolution  LIDAR  topography  data  were  collected  and  compiled  for  the 
Isaac-Verot  watershed. 


Figure  2:  An  example  of  GIS  analysis:  Soil  classifications  maps  within  the  Vermillion 

Watershed. 
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iii.Development  of  data-driven  models  for  predictions  of  complex  discharge-stage 

relations  in  low-gradient  streams. 

The  field  measurements  showed  that  development  of  stage-discharge  relationships 
for  coastal  low-gradient  streams  is  a  challenging  task.  Such  relationships  are  highly 
nonlinear,  non-unique,  and  exhibit  multiple  loops.  The  investigators  efforts  to  use 
conventional  parametric  regression  methods  to  model  these  relationships  were 
unsuccessful.  Therefore,  attempts  were  made  to  take  advantage  of  two  data-driven 
computation-intensive  modeling  techniques,  artificial  neural  networks  and  local 
nonparametric  regression.  The  results  show  considerably  better  performance  of 
both  modeling  techniques  over  conventional  regression  techniques  (Habib  et  ah, 
2005).  Both  neural  network  and  local  regression  models  were  able  to  predict  and 
reproduce  the  stage-discharge  multiple  loops  observed  at  the  outlet  of  Isaac  Verot 
catchment. 


0  500  1000  1500  2000  2500 

Discharge  (cfs) 

Figure  3:  Relationship  between  stage  and  discharge  in  flat  watersheds  subjected  to 
tidal  influence  is  complex  and  quiet  challenging.  Data-driven  models,  such  neural 
networks,  can  reproduce  these  relationships  with  relatively  high  accuracy 

iv.  Performance  evaluation  of  GSSHA  modeling  system  compared  to  other  widely 

used  hydrologic  models. 

The  purpose  of  this  analysis  was  to  gauge  the  performance  of  GSSHA  against 
MIKE-SHE  which  represents  a  more  mature  modeling  system  that  has  been  tested 
and  proved  accurate  in  numerous  research  and  applied  studies.  Both  modeling 
systems  have  similar  formulations  and  representations  of  various  rainfall-runoff 
processes.  In  addition,  a  third  modeling  system,  HEC-HMS,  was  also  included  in 
the  assessment  exercise  due  to  its  rather  simplified  and  conceptual  formulation. 
The  three  models  were  calibrated  and  validated  for  the  Goodwin  Creek  watershed 
for  a  series  of  continuous  simulation  periods  during  1982,  1987,  1997,  1998,  and 
2002.  The  three  models  were  evaluated  based  on  their  ability  to  reproduce  times 
and  magnitudes  of  runoff  peaks  and  also  runoff  volumes. 

The  overall  comparison  shows  that  GSSHA’s  results  are  highly  comparable  with 
those  of  MIKE-SHE  and  HEC-HMS.  Although  the  different  statistics  for  the  three 
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models  have  close  values,  GSSHA’s  errors  show  different  trends  from  the  other  two 
models.  Generally,  both  MIKE-SHE  and  HEC-HMS  tend  to  overestimate  the 
observed  peaks  and  runoff  volumes,  while  GSSHA  tends  to  underestimate  them. 
As  for  the  error  in  time  to  peak,  GSSHA’s  results  are  more  comparable  to  HEC- 
HMS  than  MIKE-SHE.  Both  GSSHA  and  HEC-HMS  have  mostly  positive  time  to 
peak  error,  which  indicates  a  slower  watershed  response  in  both  models 
representations  than  the  real  watershed  response  to  the  same  given  precipitation. 
Comparatively,  MIKE-SHE  had  mostly  negative  time  to  peak  error  values.  The 
RMSE  values  for  GSSHA’s  results  are  also  comparable  with  those  of  the  other  two 
models.  It  is  obvious  that  GSSHA  has  the  highest  RMSE  value  when  the  observed 
peaks  are  small  in  magnitude.  GSSHA’s  RMSE  values  become  more  comparable 
with  those  of  the  other  two  models  when  the  peaks  have  larger  magnitude. 


Date 


Figure  4:  Observed  and  simulated  hydrograph  for  the  three  models  during 

calibration  period 


Date 

Figure  5:  Observed  and  simulated  hydrographs  for  the  three  models  during  a 

validation  period 

It  is  emphasized  that  interpretation  of  the  above  results  should  take  into 
consideration  the  fact  that  the  three  models  were  calibrated  manually.  Although 
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efforts  were  made  to  maintain  as  unified  calibration  criterion  as  possible,  subjective 
judgment  on  when  to  stop  further  calibration  of  a  certain  model  was  inevitable.  An 
automated  calibration  procedure  where  all  models  are  strictly  subjected  to  the  same 
criterion  is  one  way  to  eliminate  any  ambiguities  in  interpretation  of  the  results. 

v.  Assessment  of  sensitivity  of  runoff  predictions  to  spatial  and  temporal  rainfall 
sampling. 

The  purpose  of  this  analysis  was  to  assess  the  impact  of  the  spatial  and  temporal 
sampling  resolution  of  rainfall  on  the  predictive  ability  of  both  distributed  and 
lumped  hydrologic  models.  This  analysis  is  intended  to  address  two  issues  that  are 
important  for  both  research  and  practical  purposes:  are  distributed  physics-based 
models  more  sensitive  to  degrading  the  rainfall  sampling  resolution?  And  what  are 
the  rainfall  sampling  requirements  that  can  be  considered  sufficient  for  runoff 
prediction? 

The  impact  of  the  temporal  and  spatial  sampling  of  rainfall  on  the  perfonnance  of 
physically-based  and  conceptual  models  was  investigated.  Simulations  were 
perfonned  with  varying  the  number  of  rain  gauges  from  1  to  30  gauges.  In 
addition,  temporal  resolution  of  rainfall  sampling  was  varied  from  15  minutes  to  6 
hours.  The  study  showed  that  errors  introduced  by  coarse  sampling  scenarios  can 
be  significant.  Overall,  for  the  particular  size  of  the  study  watershed,  increasing  the 
rain  gauge  density  from  1  to  2  resulted  in  the  most  significant  improvement  for  both 
models.  Similarly,  a  temporal  sampling  frequency  beyond  1  hour  showed 
significant  deterioration  in  the  quality  of  the  runoff  prediction.  This  study  also 
showed  that  the  distributed  physically-based  model  was  more  sensitive  to  the 
rainfall  temporal  and  spatial  sampling  than  the  conceptual  semi-distributed  model. 
Such  sensitivity  was  more  pronounced  and  persistent  especially  when  the  spatial 
sampling  was  significantly  lowered.  The  combined  spatial-temporal  sampling 
experiment  showed  that  increasing  the  temporal  sampling  compensates,  at  least 
partially,  for  the  loss  of  rainfall  spatial  information.  It  also  showed  that  under  poor 
spatial  sampling  conditions,  the  gain  in  model  performance  by  increasing  the 
temporal  sampling  frequency  becomes  negligible. 


Figure  6:  Error  in  rainfall  and  runoff  due  to  rainfall  temporal  resolution. 
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Figure  7:  Error  in  rainfall  and  runoff  due  to  changes  in  spatial  resolution  of  rainfall 

vi.  Using  radar  rainfall  information  for  hydrologic  modeling  with  GSSHA. 

Recent  years  have  witnessed  significant  advances  in  development  of  operational 
radar-rainfall  products.  These  products  are  desirable  for  several  hydrologic 
applications  such  as  flood  forecasting  and  rainfall-runoff  modeling.  It  is  recognized 
that  radar-rainfall  estimates  are  associated  with  unknown  uncertainties.  The  nature 
of  these  uncertainties  and  their  impact  on  the  prediction  accuracy  of  hydrologic 
models  is  a  challenging  research  problem  that  has  been  identified  by  the  hydro¬ 
meteorological  community.  Full  investigation  of  such  a  problem  is  beyond  the 
scope  of  this  project.  Therefore,  the  project  focused  on  one  aspect  of  the  radar  error 
characteristics,  which  is  systematic  error  (bias)  and  its  impact  on  runoff  predictions. 
Radar-rainfall  products  used  in  the  study  are  the  hourly,  4x4  km2  resolution  NWS 
Multi-sensor  Precipitation  Estimator  (MPE)  radar-rainfall  estimates  over  the 
Goodwin  Creek  experimental  watershed.  The  physically  based  fully  parameter 
distributed  model  (GSSHA)  was  used  to  perform  rainfall-runoff  simulations  driven 
by  different  scenarios  of  radar  and  rain  gauge  estimates.  Analysis  of  radar 
systematic  error  showed  that  MPE  radar  estimates  can  have  large  biases  in 
estimating  storm-total  rainfall  volumes.  Such  biases  cause  significant  deterioration 
in  runoff  predictions.  Bias  adjustment  over  a  stonn  time  scale  has  lead  to  a 
significant  improvement  in  runoff  prediction  accuracy.  However,  after  removing 
(adjusting)  storm  biases,  significant  random  errors  still  exist  in  the  radar 
infonnation.  Future  work  will  focus  on  quantifying  the  spatial/temporal 
characteristics  and  dependence  structure  of  random  radar  errors  and  how  their 
implications  for  uncertainties  in  rainfall-runoff  simulations. 
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Figure  8:  The  21  -km  Goodwin  Creek  Experimental  Watershed  in  northern 
Mississippi  is  selected  for  radar-rainfall  runoff  analysis  because  of  its  dense 
coverage  with  precipitation  and  surface  runoff  measuring  stations. 


Spatial  variability  of  radar  bias 


Storm  Number 


Figure  9:  MPE  radar  estimates  were  found  to  have  significant  bias  during  most  of 
the  1 1  storms  selected  for  this  study.  Storm  bias  factors  for  each  pixel  covering  the 
watershed  were  computed  as  the  ratio  of  total  radar  to  total  gauge  accumulations. 
Bias  factors  showed  significant  variability  from  one  pixel  to  another. 
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